Black Holes in the Universe: Generalized Lemaitre-Tolman-Bondi Solutions 
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We present new exact solutions which presumably describe black holes in the background of a 
spatially flat, pressureless dark matter (DM)-, or dark matter plus dark energy (DM+DE)-, or 
quintom-dominated universe. These solutions generalize Lemaitre-Tolman-Bondi metrics. For a 
DM- or (DM+DE)-dominated universe, the area of the black hole apparent horizon (AH) decreases 
with the expansion of the universe while that of the cosmic AH increases. However, for a quintom- 
dominated universe, the black hole AH first shrinks and then expands, while the cosmic AH first 
expands and then shrinks. A (DM+DE)-dominated universe containing a black hole will evolve 
to the Schwarzschild-de Sitter solution with both AHs approaching constant size. In a quintom- 
dominated universe, the black hole and cosmic AHs will coincide at a certain time, after which the 
singularity becomes naked, violating Cosmic Censorship. 
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I. INTRODUCTION 

What is the relation between cosmic expansion and 
local physics? Do local gravitational system, i.e., stars, 
galaxies, galaxy clusters, or even black holes, expand with 
the cosmic expansion? This issue has received much at- 
tention over the years [H-[2l|. however, these discussions 
often produced contradictory results leaving the issue 
rather ambiguous (see [23l - l25| for reviews). It is gener- 
ally believed that these contradictions are due to the use 
of different or unphysical coordinates, or of unphysical 
solutions |9j. 

In order to address this issue, McVittie introduced 
his renowned solution [26] in 1933; this was intended 
as describing a point mass embedded in a Friedmann- 
Robertson- Walker (FRW) universe, but modern studies 
revealed that the McVittie metric cannot describe a point 
mass 0, H, [ll|, H3| , but is rather interpreted as describ- 
ing a black hole in a cosmological background 0, |47j • 
Einstein and Strauss [28[ introduced the Swiss-cheese 
model which is, however, unable to describe the Solar 
System and suffers from other limitations fl|. [ 2 3| . Next. 



Wm, and 



the Vaidya [29j], Thakurta, Sultana-Dyer 
Faraoni- Jacques [3l| solutions were found, and these so- 
lutions also describe black holes embedded in FRW uni- 
verses, each with some restrictions. Either the cosmic 
fluid is restricted to be a cosmological constant (yielding 
the Schwarzschild-de Sitter solution, a very special case, 
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indeed), or it must be a mixture of two perfect fluids, one 
of which is a null dust (for the Sultana-Dyer solution) , or 
it is an imperfect fluid with heat flow instead of a simple 
perfect fluid (for the Faraoni- Jacques solution) [U [53] ■ 
It is of interest to look for simpler exact solutions with a 
single perfect fluid as a source to address the issue of cos- 
mic expansion versus local physics. The purpose of the 
present paper is to derive such new solutions representing 
a black hole in a FRW universe with a cosmic fluid con- 
sisting of dark matter, or of dark matter plus dark energy 
(possibly a cosmological constant), or of quintom matter 
[32j. These types of solution are also useful, if nothing 
else as toy models, in the study of dynamical horizons 
and their thermodynamics (22), l4^H45j and of primordial 
black holes [33[ as probes of the early universe [H, [35[ . 

Regarding the relation between cosmic expansion and 
local physics, Noerdlinger and Petrosian conclude their 
work [4j by stating: "Consider the possible expansion of 
clusters or superclusters of galaxies, of mean rest-mass 
density p c , immersed in a universe containing a gas of 
particles having energy density p. It is shown that when 
p > pc, the clusters or superclusters expand with the uni- 
verse, but if p c ^> p, the expansion is reduced in the ratio 
pj Pc ■ Unfortunately, this conclusion is drawn from the 
McVittie solution, which can not describe realistically a 
point mass, a star, or a cluster of galaxies in a FRW 
background. A similar conclusion was reached by Price 
[l8| who discovered the "all or nothing" behavior (i.e., 
weakly coupled systems are comoving while strongly cou- 
pled ones resist the cosmic expansion). However, Price's 
result applies only to a de Sitter background, not to a 
general FRW universe. As is well-known, de Sitter space 
is very special: the de Sitter metric can be put in static 
form in the spacetime region between the black hole and 
the cosmological horizons, which explains why strongly 
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bound systems in this background do not expand 31]. 

In this paper, we present new solutions which we pro- 
pose to interpret as describing cosmological black holes 
embedded in FRW backgrounds, with the advantage over 
previous solutions [HI, 0, [29j - [3l| that they are sourced by 
a single perfect fluid, without heat flux. These solutions 
lead to surprising results: the black hole AH decreases 
with the expansion of the universe when the latter is 
DM- or (DM+DE)-dominated. Since the Misner-Sharp 
mass is a half of the radius of the AH, also the black hole 
mass decreases with the cosmic expansion. As the most 
strongly bound gravitational system, it is intuitive that 
the size of a black hole should increase with the expan- 
sion of the universe due to the swallowing of surround- 
ing cosmic matter. However, this picture is incorrect in 
general as the cosmic expansion "wins" over the local 
gravitational attraction of the black hole. The universe 
is always expanding after the Big Bang and the energy 
density of the cosmic fluid decreases with cosmic time. 
The seemingly bizarre behavior of the Misner-Sharp mass 
derives from the fact that it is really a mass sum (i.e., 
the mass of the background fluid is included) , and coin- 
cides with the Hawking-Hayward quasi- local mass [2lj |. 
We note that this is a purely classical effect which makes 
the black hole mass decrease and has nothing to do with 
Hawking radiation [36]. 

In a quintom-dominated universe, during the mat- 
ter era the black hole mass decreases, while during the 
phantom-dominated epoch it increases. The physical rea- 
son is that the cosmic density first decreases and then 
increases. Since the Misner-Sharp (MS) mass is the sum 
of the black hole mass and of the background mass, the 
black hole mass first decreases and then increases. One 
interesting result is that the black hole singularity will 
become naked before the Big Rip occurs, which violates 
the Cosmic Censorship Conjecture 37). So if the latter 



coordinates as 



is correct, phantom matter may not exist in nature. 

In this paper, we use units in which the speed of light 
c and Newton's constant G assume the value unity, and 
the signature +2 for the metric. The paper is organized 
as follows: first we show that several known solutions 
can not describe a black hole embedded in a matter- 
dominated universe. Second, we present new solutions 
describing a black hole in a matter-dominated universe, 
a dark matter plus cosmological constant-dominated uni- 
verse, and a quintom-dominated universe. The last sec- 
tion contains the conclusions. 



II. THE MCVITTIE SOLUTION 

A natural starting point to investigate the relation be- 
tween cosmic expansion and local physics is the 1933 
McVittie solution [26J, which can be written in isotropic 
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where the constant M reduces to the physical black hole 
mass when a(t) =const. (this can be seen by rescaling 
the radial coordinate r — > r = ar), and the scale fac- 
tor a(t) is an arbitrary function of the cosmic time t. 
When M = 0, the solution reduces to the spatially fiat 
FRW metric, and at a first glance, it can be understood 
as representing a Schwarzschild black hole embedded in 
a spatially flat FRW universe |4(. However, it is now 
known that, with the exception of the Schwarzschild-dc 
Sitter solution, matters are more complicated because 
the McVittie metric is singular on the 2-sphere r = 2Mo<x 
(which reduces to the Schwarzschild horizon if a =const.) 
Q and this singularity is spacelike [ll|. It was claimed 
in the past that the McVittie metric describes a point 
mass located at r = and embedded in a FRW universe. 
However, this point mass is, in general, surrounded by 
the singularity at r = Mo / (2a) . This singularity was 
studied in Ref. [6]. Nolan Q showed that it is a weak 
singularity in the sense that an object falling onto the 
r = Mq/{2o) surface is not shrunk to zero volume, and 
therefore the energy density of the surrounding fluid is 
finite. However, the pressure diverges at r — Mq/ (2a) to- 
gether with the Ricci scalar R [(| [27j . While there is little 
doubt that the McVittie metric represents some kind of 
strongly gravitating central object, its physical interpre- 
tation is not completely clear and is still a subject of 
debate jllllSi, S3- 

In view of this situation we show that, with the excep- 
tion of the Schwarzschild-de Sitter solution, the McVit- 
tie solution is unable to describe a black hole or a point 
mass in a FRW universe when the equation of state of 
the cosmic perfect fluid is p — w(t)p. The proof of this 
statement proceeds as follows. The McVittie solution 
was obtained by forbidding explicitly the accretion of 
cosmic fluid onto the central object. This requirement 
corresponds to Gj = 0, which in turn implies that the 
stress-energy tensor component Tq 1 = and there is no 
radial flow of energy. Let us assume that the energy- 
momentum tensor that sources the McVittie metric is 
that of a perfect fluid, — (p + p) U^JJ V + pg^ v , where 
is the 4-velocity of the comoving observer, and p and 
p are the density and pressure of the fluid, respectively. 

By substituting the metric and the energy-momentum 
tensor into the Einstein equations — 8ttT^ u , the re- 
sulting Einstein equations can be arranged as 



3H 2 



8np , 



(2) 



^+ff 2 = -8np-^(H-*Hi-4np) , (3) 
a ar V 2 / 
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where H = a/a is the Hubble parameter and an overdot 
denotes differentiation with respect to the cosmic time t. 
We recognize eqs. j2j and (j3|) as the Friedmann and the 
acceleration equations, respectively. 

We have p = and p — p/3 for the matter-dominated 
and the radiation-dominated eras, respectively, while 
p = w{t)p with w < —1/3 for a dark energy-dominated 
universe. For all these cases, the equation of state as- 
sumes the form p = w(t)p. Eq. ([2]) tells us that p depends 
only on t and, since p = w(t)p, the same is true for p. In 
eq. ((3j, the expressions — + H 2 and —inp depend only 

on t, while — LH - | H 2 - 4npj depends on both t 

and r, therefore eq. (|3]) is satisfied if and only if 



of the McVittie solution 



3H 2 = 8np, — + H 2 + 8np = 
a 



H — -H — 4-rrp = . 



(4) 



These admit the unique solution 

H = const. , p = —p = const. , a(t) = e Ht , (5) 

which is the Schwarzschild-de Sitter spacetime. There- 
fore, the McVittie solution can not describe a black hole 
embedded in a matter- dominated, radiation- dominated, 
or dark energy- dominated universe with perfect fluid 
equation of state p — w(t)p. In one word, given the equa- 
tion of state p — w{t)p, the McVittie solution is nothing 
but the Schwarzschild-de Sitter solution. 

Actually, Kaloper et al. 46| have argued that in the 
presence of a cosmological constant, McVittie's solution 
does in fact represent a black hole in an expanding uni- 
verse. Here we provide independent support for this 
statement. 

Recently, assuming a specific function H(t) in the 
McVittie solution, Lake and Abdelqader [47| performed 
a detailed study of the particular McVittie solution and 
corroborated the conclusion of Kaloper et al. Moreover, 
using a tetrad-based method for solving Einstein's field 
equations, Nandra et al. [48| obtained solutions describ- 
ing objects embedded in FRW universes which provide 
a new perspective on the McVittie metric. The study 
of the global structure of the McVittie solution in |47| 
revealed unexpected features which make this solution 
much more complicated than it would seem at a first 
sight. Something similar probably happens for the new 
solutions presented here, although we will leave the de- 
tailed investigation of radial null geodesies and of the 
global structure to future work. 



III. THE FARAONI- JACQUES SOLUTION 

Because the McVittie solution can not describe a black 
hole in a FRW universe, we turn our attention to the 
Faraoni- Jacques solution [3lj |. which is a generalization 
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(6) 



in which the function M(t) replaces the constant Mq 
of the McVittie metric. This corresponds to lifting the 
McVittie non-accretion restriction Gq = 0. This solution 
seems to describe a cosmological black hole embedded 
in a FRW background and presents advantages over the 
previous Thakurta and Sultana-Dyer solutions [HI, [3(| , 



ds 2 = - 



1 



T dt 2 



(dr 2 



2 dfl 2 



(7) 



in the sense that both the energy density and the pressure 
are finite near the black hole horizon, and the energy den- 
sity is positive-definite [3l[ . Now let us examine whether 
it can describe a black hole embedded in FRW universe 
sourced by a perfect fluid. 

We assume that the cosmological matter is described 
by the single perfect fluid energy-momentum tensor 



Tfu, = (p + p) U^Uv + pg^y 



(8) 



where is the 4-velocity of the fluid. The only non- 
vanishing components of the Einstein tensor are 



g\ = g\ = G\ 7^ . 



(9) 



If a radial energy flow is allowed, the fluid four- velocity 
is 



and the normalization U^U^ 



IP* = (C/°, L/ x ,0, 0) , 

-1 gives 



a 2 v 2 g m \ 1 



M{t) 
2a (t) r 



,v,0,0 



(10) 



,(H) 



where v is the proper 3-velocity of the fluid. From Gn = 
G22 we conclude that 



(p + p)a 2 v 2 1 + 



M(t) 



= 



2a (t)r, 

We have two types of solutions characterized by 

v = , M = M 
(which is the McVittie solution), or 
P = ~P, 



(12) 

(13) 
(14) 
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which gives the Schwarzschild-de Sitter solution. There- 
fore, the Faraoni- Jacques solution (and also the Thakurta 
and Sultana-Dyer solutions) can not describe a black hole 
embedded in a matter-dominated universe. In general, 
they can not describe a black hole embedded in a single 
perfect fluid- dominated universe. However, the Faraoni- 
Jacques solution can describe a black hole embedded in 
FRW Universe provided that the source is given by an 
imperfect fluid with a heat flux: 



(p + p) U^U V + pg^y + q ll U v + q^Ufj, 



(15) 



where = [^-g m , 0,0,0 and q» = (0,g,0,0), where 



q is the heat flux density. This has been shown by us in 



IV. THE VAIDYA SOLUTION 

The Vaidya metric is obtained by conformally trans- 
forming the Minkowski one: 

ds 2 = a 2 (rj) (-dn 2 + dr 2 + r 2 d6 2 + r 2 sin 2 



— (dr/ + dr) 2 



(16) 



where Mq is the black hole mass and 77 is the conformal 
time. When Mq = and a = n 2 , this metric reduces 
to that for a matter-dominated universe. On the other 
hand, when a =const., it reduces to the Schwarzschild 
metric. At a first glance, one is tricked again into be- 
lieving that this metric may describe a black hole in a 
matter-dominated universe but this is incorrect, as shown 
below. 

The energy-momentum tensor of dust or dark matter 

is 



PdUJJv 



(17) 



Due to spherical symmetry, the four-velocity is writ- 
ten as eq. (p~0|) . We have two vanishing components of 
T av , i.e., T£ and T 3 3 , hence G\ = G\ = 0. Given the 



metric (16]), G\ 



2a 



G3 is calculated as [27| . 



/ 2md 4ma z 
\ ra 5 ra 6 



= 



(18) 



The terms in the first bracket depend only on time, while 
those in the second bracket are r-dependent, so eq. (fT8|) 
is satisfied if and only if 



2d 



a a a* a a D 



(19) 



Then we necessarily have the trivial solution a =const., 
which is just the Schwarzschild metric. This shows that 
the Vaidya solution can not describe a black hole embed- 
ded in a matter- dominated universe, cither. 

Thus far, we have shown that several known exact solu- 
tions can not describe a black hole in a matter-dominated 
universe. In the next section, we look for such a solution. 



V. BLACK HOLE IN A MATTER UNIVERSE 

A spherically symmetric inhomogeneous gravita- 
tional field produced by dark matter with energy- 
momentum (|17l) is described by the Lemaitre-Tolman- 
Bondi (LTB) metric [38| which, in the notation of [39j j . 
is given by 



ds 2 



-dt 2 



R' 



1 + / 



dr 2 + R 2 dn 2 



(20) 



Here / is an arbitrary function of the comoving coordi- 
nate r satisfying / > 1, R(t,r) is the physical radius at 
time t and coordinate radius r, while a prime represents 
differentiation with respect to r. 

With the energy momentum tensor (Tl"T|) and the met- 
ric (f2"0"|) . the Einstein equations read 



3(i? 2 -/ 
"R 2 

F' 
R 2 R' 



8ttF 



3pd 



(21) 



(22) 



where an overdot and a prime represent partial differen- 
tiation with respect to t and r, respectively, and F(r) is 
an arbitrary function of r. If we substitute 



/ = -for 2 , 
into eqs. ([2"Tjl and 



F 



■ Por- 



R = a(t)r 



the latter become 



3 (a 2 + f ) _ 8n Po 



Pd 



Po 

«3 



(23) 



(24) 



with /o and po constants. We recognize eq. (|24|) as the 
Fricdmann equation with a dust source, hence /o rep- 
resents the spatial curvature and po the density of the 
universe at time do = 1. The observed universe is spa- 
tially flat with good accuracy (probably due to inflation) 
and we can neglect the curvature term in what follows 
by setting / = 0. 

For an arbitrary function F(r) and / = 0, the solution 
of eq. (|2ip was given by Tolman 40] and Oppenheimer 
and Snyder (4lj in their pioneering investigations of the 
gravitational collapse of dust as 



R = 



, / x 3 /8tt „ 



2/3 



(25) 



where the integration "constant" h is an arbitrary func- 
tion of r. Substitution of eq. ([25]) into ([20)) yields 



ds 2 



-dt 2 



R- [ ' !h 



6?r td VF 



R 2 dQ 2 . (26) 



There is actually less freedom in eq. (|25p than is apparent 
from the two arbitrary functions F and h since F can be 
taken as a new radial variable. It is convenient to choose 



,3/2 



(27) 
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and 



R = 



2/3 



(28) 



Eq. (|2"8")l is the general solution representing a spherically 
symmetric gravitational held with dust as a source. As 
shown above, this can describe a dust-dominated uni- 
verse. On the other hand, if the dust disappears, it de- 
scribes a vacuum, spherically symmetric, gravitational 
field and reduces to the Schwarzschild metric due to 
the Jebsen-Birkhoff theorem. By combining these con- 
siderations, one guesses that eq. (f28f may describe a 
Schwarzschild black hole immersed in a dust-dominated 
universe. In the following, this claim is examined. 



A. Minkowski space 

By solving the Einstein equations with p e i = 0, one 
finds the vacuum solution 



F = F 



(29) 



where Fo is an integration constant. If Fo = 0, we obtain 
R = r and we recover the Minkowski solution. 



i.e., the Schwarzschild solution. The physical meaning 
of Fo is derived from Fo = ^/2rn and the Schwarzschild 
solution is given by 



R = 



'2m t 



2/3 



(35) 



C. Dust-dominated universe 

Eqs. (f22|) and (f28|) show that the homogeneous and 
isotropic, spatially flat, dust-dominated FRW solution is 
recovered if and only if F = pof 3 . The dust universe 
solution is then given by 

, -, 2/3 



R 



(36) 



The scale factor and the density are 

2/3 



and 



\J 67T/90 t + 1 



Po 



(37) 



(38) 



respectively, and it must be t > to, with to 
representing the Big Bang. 



B. The Schwarzschild solution 

If Fo ^ 0, we recover the Schwarzschild solution. It is 
convenient to set F = with Fo a positive constant. 
Then, the Schwarzschild solution is given by 



R = 



,3/2 



-i 2/3 



Fnt 



(30) 



Fo is determined by the mass m of the black hole. 
In order to show this, let us rewrite the metric in 
Schwarzschild coordinates, setting 



x(t,r) 



r 3/2 + - F t 



2/3 



(31) 



Then, eq. (|20p becomes 
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2F 



dtdx + x^dQ, 2 . (32) 



Using the new time variable 



T = t - t + 2F V^ + ^o 2 ln 



Vi-Fo 
Vi + Fo 



(33) 



one obtains 



ds z 



' P2 - 

l-t±\dT 2 



1 - — ) dx 2 + x 2 dn 2 



(34) 



D. Black hole in a dust-dominated universe 

By combining the considerations above, one expects 
the metric for a black hole immersed in a dust universe 
to be given by 

-i 2/3 



R 



,3/2 



'2m t - 



87T P0 r 3/2 ^ 



(39) 



If m — 0, this line element describes a dust universe, 
while if po = 0, it reduces to a Schwarzschild solution; 
(|39l) should describe a black hole immersed in dust uni- 
verse. It is clear, by comparison of eqs. (f3"9"]) and (|2"5j) . 
that this is a solution of the full Einstein equations. The 
substitution of eq. (|3"9")l into eqs. (|2"Tj) and (|2"2"|) gives the 
energy density 



Pd 



po"(3V7^+2V3w 3 / 2 ) 



2 % /3^( % /6Vp "i+ l)i?3/2 



(40) 



The physical radius R cannot be negative, hence we re- 
quire that r > and t > to = — ^/Girpo- Then, the density 
is always positive. If m = 0, we have pd = po/R 3 , i.e., 
dust or dark matter. Since R = represents the Big 
Bang, the latter occurs at the time 



t 



-V2mr^ 3/2 
2 



\/6ttpo 



(41) 



and the energy density is positive everywhere on the 
spacetime manifold. 
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VI. EVOLUTION OF THE APPARENT 
HORIZONS 



where J (t, x) is an integrating factor which always exists 
and solves the partial differential equation 



In the Schwarzschild-de Sitter spacetime, there exist 
a black hole apparent horizon (AH) and a cosmic AH. 
The metric (|20[) describes a spherically symmetric and 
dynamical black hole in a FRW background more gen- 
eral than de Sitter. In this case, the event horizon may 
not be well-defined, but the AH always exists. The AH 
is a marginally trapped surface with vanishing expansion 
and has been argued to be a causal horizon for a dynam- 
ical spacetime. The AH is associated with the Hawking 
temperature, gravitational entropy and other thermody- 
namical aspects |42h45| | . The first law of thermodynamics 
for the AH has been derived not only in general relativity 
but also in several other theories of gravity including the 
Lovelock, nonlinear, scalar-tensor, and braneworld theo- 
ries [49-55]. In view of this point, in order to investigate 
the evolution of the black hole mass, we calculate the 
energy contained inside the AH which has been used in 
black hole thermodynamics in relation with energy flows 
through the AH [49|, [56( . Let us proceed to calculate the 
radius of the black hole AH. 

For a spherically symmetric spacetime with line ele- 



h^dx^dx" 



ment ds 2 

Sharp mass is Mms 
AH, it is W v 
mass inside the AH is simply 



dfl , the generalized Misner- 
= x (1 - W v x tli x, v ) /2 [13. At the 
and the generalized Misner-Sharp 



M MS = xah/2 



(42) 



where xah is the radius of the AH. In order to know the 
evolution of the black hole MS mass, we need to know 
the AH radius. 

We rewrite eq. ([2"0"| in Schwarzschild coordinates using 
the new spatial coordinate 



x 



2/3 



(43) 



instead of r, in terms of which eq. (|20)) is rewritten as 

2" 



ds 2 = - 




dt 2 + dx 2 



_ JL (^/2m~ + ^pr 3 ^ dtdx + x 2 dQ 2 . (44) 

The coordinate system (t, x, 9, 4>) is not orthogonal but 
the cross-term dtdx can be eliminated by introducing the 
new time T defined by 



dT : 




dxJ- 1 = d t I -J- 1 



1 




(46) 



which guarantees that dT is an exact differential. The 
metric (l4"4")l is cast in the Schwarzschild form 



ds 2 




J 2 dT 2 



dx 2 



(47) 



and the equation of the AHs is 
f 




W r 3/2 ] _ q 



(48) 

Substituting eq. (|4"3")l into eq. (|48p. we obtain 

2y/6poTTx 3/2 - y/x (3 + 3^67r p tj + Z\p2m = . (49) 

This cubic equation for yfx has, in general, only two pos- 
itive roots which represent the black hole AH and the 
cosmic AH, respectively. In fact, if po = 0, we obtain 
from eq. (|49]) 



x = 2m , 



(50) 



the AH of the Schwarzschild black hole. On the other 
hand, if m = 0, it is 



X= 2 



1 



671-po 



+ t 



(51) 



which is the AH of a dust-dominated universe. The black 
hole and the cosmic AHs are plotted in Fig. [1] As is 
clear from this figure, the radius of the black hole AH 
decreases while that of the cosmic AH increases as the 
universe evolves. There was an early time at which the 
two horizons coincided, and before which both horizons 
were absent and the singularity was naked. Since the MS 
mass is proportional to the radius of the black hole AH, 
the MS mass of a black hole decreases with the evolution 
of the universe. The reason for this behavior can be un- 
derstood as follows: with the expansion of the universe, 
the cosmic density is decreasing, which is equivalent to a 
decreasing of black hole AH, so the black hole MS mass 
decreases. 
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FIG. 1: The size of the black hole AH (lower curve) in a 
dust-dominated universe decreases with the cosmic expansion 
while that of the cosmic AH (upper curve) increases. The MS 
mass of the black hole tends to zero in the future. There 
was a time at which the two horizons coincided and before 
which they were absent and the singularity was naked. The 
plots correspond to the parameter values m = 1, po = 0.05 
in Plank units. In the International Units System, they are 
m = 10 24 M Q and p = 0.05 ■ 10 -123 p P (p P is the Plank 
density) . 




0.2 0.4 t 0-6 0.8 1 



FIG. 2: The AH of an astronomical black hole in a dust- 
dominated universe decreases with the cosmic expansion. The 
plots correspond to the parameter value m — 1O 7 M0. The 
units of x and t are 2m (the corresponding Schwarzschild 
radius) and Hq , respectively. 



We note that Fig. [T] is plotted in Plank units G = c = 
h = 1 . Translating into International units and taking 
into account the present cosmic density po = 10 -123 pp 
(where pp is the Plank density), we have a very super- 
massive black hole m = 10 24 Mq for Fig. [T] The corre- 
sponding units of x and t are 10 5 Mpc and the Hubble 
time Hq 1 , respectively. To our knowledge, such super- 
massive black holes may not exist in the universe, so in 
Fig. [2] we also plot the evolution of the AH for an astro- 
nomical black hole with mass m = 10 7 Mq. We have 
shifted the Big Bang time to t = 0. 

The naked singularity is present before the appearance 
of the two apparent horizons and cannot arise from reg- 
ular initial data. 



VII. BLACK HOLE IN A ACDM UNIVERSE 

Astronomical observations show that t he p resent-day 
universe is dominated by dark energy [HD, [5!| , hence we 
seek exact solutions describing a black hole immersed in 
a mixture of dark matter and dark energy. The corre- 
sponding stress-energy tensor is 

= pdU^Uv + pA9^v ■ (52) 

The first term on the right hand side describes dark mat- 
ter and the second term describes dark energy. p\ is a 
positive constant 1 We assume the line element to be of 
the form 

ds 2 = -dt 2 + e^dr 2 + e^dn 2 , (53) 

where Q and w are functions of t and r. In comoving 
coordinates, the four-velocity is U 1 * — (1,0,0,0) and the 
Einstein equations are 

G° = 87r( Pd + p A ) , G? = 0, 

G\ = 87rp A , G\ = 871-pA . (54) 

From G\ = one obtains 

e Q = e"w' 2 /4 . (55) 

Note that there is no accretion onto the central object, 
as in the McVittie metric. Following Oppenheimer and 
Snyder, we have set the integration "constant" (f(R) in 
41]) to unity. Substituting eq. (f5"5"T) into eq. (|5"5|). we 
obtain 

3 

u)' + - ljlo' = . (56) 

Integration with respect to r yields 

w + ^ih 2 = -8%K(t) , (57) 

where K is an integration "constant" (an arbitrary func- 
tion of t). This equation is consistent with G\ — 8np\ if 
and only if 

K (t) = -p A . (58) 
The Einstein equations then reduce to 

G\ = lv + -d} 2 = 8np\ , 

Gl = ^ + h 2 = 8ir(p d + PA ) . (59) 

id 4 



There is evidence in [46l | and [47| that the presence of a cos- 
mological constant makes a significant difference in the McVittie 
solution, and we expect that this is the case also for our gener- 
alized Lcmaitre-Tolman-Bondi solutions. 
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Using k = y/6np\, the Einstein equations have the 
solution 



e" = (Se kt + Pe- kt ) 4/3 



(60) 



where S and P are arbitrary functions of r. Substitu- 
tion of eqs. (l60|) into eq. (l53|) and the use of S as a new 
spatial coordinate leave only one degree of freedom. It is 
convenient to choose 



S = r 3 ' 2 , 



(61) 



which has the advantage that if P = and p\ = 0, the 
line clement (53) reduces to the Minkowski one. 



A. The Schwarzschild-de Sitter solution 



Comparing eq. (|67[) with the well-known form of the 
Schwarzschild-de Sitter solution 



ds 2 = f 1 - *2 - ^x 2 I dT 2 + 



x 3 

one deduces the physical meaning of Sq, 

9m 



(68) 



(69) 



Then, the Schwarzschild-dc Sitter solution in comoving 
coordinates is given by 



3/2 kt _ ^ -3/2 -kt 
8fc 2 



4/3 



e<V 2 /4 . 



(70) 



The Einstein equations with pd = yield 

P = S r~ 3/2 , (62) 

where So is an integration constant. The Schwarzschild- 
de Sitter solution is given by 



= (r 3 ' 2 e kt + Sbr- 3 / 2 e - fet ) 



4/3 



e"w' 2 /4 



(63) 



The physical meaning of So can be understood after 
rewriting eq. (|63[) in Schwarzschild coordinates; setting 



x = (r 3 ' 2 e kt + 5 r- 3 / 2 e- fct ) 
the line clement (1551) becomes 



2/3 



(64) 



ds 2 = - 



4k 2 (x 3 - 4S ) 
9x 



dt 2 + dx 2 



4k 



3y/x 

Introducing the new time 



y/x 3 - 45 dtdx + x 2 dtl 2 



(65) 



T = t 



4k 2 (x 3 - 4S ) 
9x 



2k 



the metric ([631) simplifies to 



\/ x 3 — 4Sq dx 



(66) 



ds 2 



16^0 4fc 2 2 \ 2 2 

1 H ar | dT + a; as I 

9x 9 



16fc 2 So 4fc 2 



9x 
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(67) 



B. ACDM universe 



Eqs. (|60p and (|6ip show that we recover the homoge- 
nous and isotropic, spatially flat universe with a mixture 
of dust matter and dark energy if and only if we set 
P = —Tzrjr 3 / 2 . This solution is obtained for 

16k. 



r 3/2 e kt 



e UJ uj' 2 /4 



3po 
16k 2 



r 3/2 e -kt 



4/3 



(71) 



The density of dust/dark matter and the scale factor are 
given by 



Po 
„3' 



16fc 2 



2/3 



(72) 



where t > ^ In , and where t = In represents 
the Big Bang singularity. 



C. Black hole in a ACDM universe 

The inspection of eqs. (|T0"|) and (fTTj) suggests that the 
solution for a black hole in a ACDM universe should be 
given by 



9m 



£ u = ( r3 /2 e kt _ ™ -3/2,,-** _ JVO r 3/2 e -fct 

8fc 2 16/c 2 



e° = e^' 2 /4 



4/3 



(73) 



It is clear, by comparison of eqs. (f73j) and (|60|) . that this 
is indeed a solution. The substitution of eq. (|73[) into the 
Einstein equations gives the density of dust 

p Q (32r 3 e kt + 6me~ kt - p r 3 e- kt ) 2 , x 

Pd = — ■ (74) 

i? 3 (32r 3 e fet — 6me~ kt — por 3 e~ kt ) 
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This quantity is positive on the entire spacetime mani- 
fold. If m = 0, we have pd = po/R 3 , the usual scaling. 



VIII. BLACK HOLE EVOLUTION IN A ACDM 
UNIVERSE 

Here we investigate the evolution of a black hole AH 
in a ACDM universe. For this purpose, we rewrite the 
metric in Schwarzschild coordinates using 



_ p -2kt/Z r I e 2kt 



9m 3po 



16A; 1 



2/3 



, (75) 



eq. ([53]) becomes 



ds 2 



x 9 3x 



, ,'2m 4fc 2 2 po 3 \ 1/2 , , 

-2| 1 x 2 + —r 3 I dtdx 

x 9 3cc 



+dx 2 + x 2 dfl 2 . (76) 
Introducing the new time variable T with 



■ i , 2m 4fc 2 o Po r> 

dT = J~ x dt + J ( 1 x 2 - — r 3 

x 9 3x 



x 9 2>x 



(77) 



where J(t, x) is an integrating factor, we obtain a black 
hole solution in a ACDM universe. 



ds 2 





2m 


4fc 2 2 

- 

9 




a; 


3x 




2m 

:r 


4fc 2 2 


PO 3 

to' 



dx 2 



The equation of the AHs is 



so 



(78) 



(79) 



(80) 



The substitution of eq. (|80j) into eq. ([To]) yields the equa- 
tion of the AHs, using which the AHs evolution is plotted 
in Fig. 3. The black hole AH shrinks while the cosmic AH 
expands with the expansion of the universe. There was a 
time in the past at which the two horizons coincided and 
before which no horizon existed and the singularity was 
naked. Late in the history of the universe, both black 
hole and cosmic AHs approach constant values. 



2m Ak 2 2 po 3 
3oo = 1 7T-x ~ ^~r° = 

a; 9 Sx 

'3a; 6m 4fc 2 x 3 \ 1 ^ 3 
Po Pa 3p 



X 

2.31 z 

2.15 z 

2.00 

1.85 

1.69 

1.54 

1.38 

1.23: 




-6.6 -3.3 t 3.3 



FIG. 3: In a ACDM universe, the black hole AH (lower 
curve) shrinks while the cosmic AH (upper curve) expands 
with the cosmic expansion. There was a time at which the 
two horizons coincided and before which no horizon existed 
and the singularity was naked. Late in the history of the 
universe, both black hole and cosmic AHs approach con- 
stant size. The plots correspond to the parameter values 
m = 2.2 • 10 22 A/ O ,p = 0.27pc,PA = 0.73pc- pc is the 
current cosmic density. The unit of t is the present Hubble 
time Hg 1 . 



IX. BLACK HOLE IN A QUINTOM UNIVERSE 

The first year WMAP data combined with the 2dF 
galaxy survey and the supernova la data favor the phan- 
tom energy equation of state of the cosmic fluid w < 
— 1 over the cosmological constant (w — —1) and the 
quintessence field (w > —1) [gFJ, [6lj]. The data seem to 
slightly favor an evolving dark energy with w < — 1 at 
the present epoch and w > — 1 in the near past [62] . The 
dark energy candidate with evolving w is named quin- 
tom [32|]. In this section, we look for the solution for a 
black hole in a quintom-dominated universe. We take the 
energy-momentum tensor for the quintom fluid to be of 
the form 



T M „ = (p + p) U^Uy + pg^ v , 



(81) 



and the line element in the form of eq. (|53[) . In comoving 
coordinates, the four-velocity is = (1, 0, 0, 0) and the 
Einstein equations are 

G° = 87175, G? = 0, 

G\ = -8tt P , G\ = -8tt P . (82) 

From Gi — (again, there is no accretion), we obtain 

e° = e"V 2 /4 . (83) 

Again, we set the integration "constant" to unity and, 
inserting eq. (|53"j) into eqs. (|52")l . we obtain 



and 



Ul + — LULO = U 



u + jilj 2 = -8irK (t) , 



(84) 



(85) 
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where K(t) is an integration "constant" . This equation is 
consistent with G\ = —87175 only if the pressure satisfies 



p = K (t) , 



(86) 



i.e., it is spatially homogenous, which can be seen as the 
consequence of three facts. First, we have set f(R) = 1 
for a spatially flat background. Second, we use comov- 
ing coordinates and, third and most important, we have 
taken the source to be a single perfect fluid. Therefore, 
the comoving observer will see a homogenous pressure. 
The Einstein equations simplify to 



u+-uj 2 = -8tt P (t) , 



— + -uj 2 = 8np(t,r) 

U! 4 



(87) 



Now we have three functions and p but only two 
equations and the system is not closed. For simplicity, 
we assign the pressure in the form 



Po 



(to - ty 



(88) 



where po is a positive constant to keep the pressure al- 
ways negative, and to is a positive constant identifying 
the Big Rip singularity. The general solution of eqs. (|87|) 
is then given by 



P(r) (t -t)^ 1 - VI ™^) 
+S(r) (t -t)^ 1+ ^™^)] 4/i 



(89) 



where S and P are arbitrary functions of r. It is conve- 
nient to choose 



p = r 3/2 



(90) 



then the candidate solution for a black hole in a quintom 
universe is 



.3/2 ( io _£)MWl+247rpo) 



y^V/2 ( tQ _ ^(l + Vl+Mirpo) 



4/3 



(91) 



If m = 0, we get a quintom-dominated cosmology while, 
if po = 0, we get a candidate solution for a black hole in 
a matter-dominated universe. Substituting eq. (|91[) into 
eqs. (|87l) . we find that the energy density is positive and 
inhomogeneous . 



X. EVOLUTION OF A BLACK HOLE IN A 
QUINTOM UNIVERSE 

To investigate the evolution of the black hole AH, we 
rewrite the metric in Schwarzschild coordinates using 

„3/2 (t _t)K 1 -^ T +24^) 
~V2^(t Q -t)i( 1+VTTm ^) 

.3 /^ r 3/ 3(to _ t) i(i+vrPH5w) m] 
2 V 3 



and define k = + 2AnpQ. Eq. (|53|) becomes 

ds 2 = - (1 - £f 2 ) dt 2 + dx 2 - 2£fdtdx + x 2 dfl 2 , 



Introducing the time coordinate 



dT = ,r 



dt- 



1-J^ 2 



dx 



(94) 



where J is an integrating factor, the metric becomes 

1 



ds 2 = -(l-£f 2 )j 2 dt 2 



dx 2 + x 2 dn 2 (95) 



1-J*f 2 



and the AHs are identified by 

1 - Jz* 2 = . 



(96) 



Eq. <(92]) then yields 



,,3/2 



1 

Wo 



2x 3 / 2 



t) 



|(i+fc) 



i(i-fc) 



- \767rpo (io - 2 



i(i+fe) 



.(97) 



Substituting eq. (|9"T|) into eq. an explicit equation 

for the AHs in terms of t and x is obtained. The evolu- 
tion of the AHs is reported in Fig. 21 In the quintom- 
dominated universe, the size of the black hole AH first 
decreases in the matter era and then increases during the 
phantom-dominated epoch (63[. The size of the cosmic 
AH, instead, first increases and then decreases. There 
exists an instant of time in the past at which the two 
horizons coincided. Both horizons disappear before the 
Big Rip, leaving behind a naked singularity and violating 
Cosmic Censorship [37| • This result is consistent with our 
previous discussion in (2lj . The naked singularity does 
not arise from regular Cauchy data, in fact it was al- 
ready present before the appearance of the two apparent 
horizons. 

The reason for the increase of the black hole mass in the 
phantom-dominated epoch can be understood as follows: 
the density of the universe is increasing, which makes the 
black hole AH also increase. Hence, the black hole mass 
(one half of the AH radius) increases. On the contrary, 
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7.5: 




-0.4 -0.2 0.2 0.4 0.6 
t 

FIG. 4: In a quintom-dominated universe, the black hole AH 
(lower curve) first shrinks and then expands, while the cosmic 
AH (upper curve) first expands and then shrinks. There exist 
two instants of time, one in the past and one in the future, at 
which the two horizons coincide. Before or after these critical 
instants, no AH exists and a naked singularity is present. 
The plots correspond to the parameter values m = l,po = 
0.0002, po = 0.001, and to = 1. 

in a matter-dominated universe, the cosmic density is 
decreasing, so the black hole AH and the black hole mass 
decrease. In a cosmological constant-dominated universe, 
the cosmic density is constant and the black hole mass 
also stays constant. 

XI. DISCUSSION AND CONCLUSIONS 

We have shown that several well-known and less- known 
solutions of the Einstein equations can not describe a 
black hole embedded in a matter-dominated universe 
sourced by a single perfect fluid. This includes the 
McVittie, Thakurta, Sultana-Dyer, Vaidya, and Faraoni- 
Jacques solutions. Motivated by this fact, we have con- 
structed exact solutions of the Einstein equations pur- 
porting to represent a black hole embedded in a back- 
ground universe which are as simple as possible, begin- 
ning with solutions with a single perfect fluid. These new 
solutions generalize the Lemaitre-Tolman-Bondi metrics 
(which are restricted to dust universes) and do not ex- 
hibit accretion of the surrounding cosmic fluid onto the 
central black hole. We have found metrics which pre- 
sumably describe a black hole immersed in a matter- 
dominated, a matter plus dark energy-dominated uni- 
verse, and a quintom-dominated universe, respectively. 

The AH and the Misner-Sharp mass are related to the 
thermodynamics of dynamical black holes and FRW uni- 
verses [42M45L |49H56| . and, therefore, we investigate the 
evolution of the black hole and cosmic AHs, as well as the 
Misner-Sharp mass. We find that the black hole mass de- 
creases in a matter- or matter plus cosmological constant- 
dominated universe. In a quintom-dominated universe, 
the black hole mass decreases in the matter era, while it 
increases in the phantom-dominated epoch. The physi- 
cal reason is that the cosmic density first decreases and 



then increases, which makes the radius of the black hole 
AH first decreasing and then increasing. Then, also the 
black hole mass (one half of the AH radius) first decreases 
and then increases. An interesting result is that the AHs 
will disappear and the singularity will become naked be- 
fore the Big Rip is reached, in violation of the Cosmic 
Censorship Conjecture [131 • If the latter is correct, then 
phantom matter may not be allowed to exist in nature. 
This statement, however, must be taken with a grain of 
salt: in fact, it is based on a particular solution of the field 
equations which may be very special. While the result 
is still interesting because very few exact solutions are 
known to represent the physical situation under study, 
the solution proposed may still be too special, or fine- 
tuned, to draw general conclusions. The phenomenology 
reported here, however, matches that obtained with very 
different classes of solutions (which are indeed accreting 
cosmic fluid) in plj . 



Another rather surprising result is that the black hole 
mass decreases in a matter-dominated universe. As the 
most strongly bound gravitational system, it seems in- 
tuitive that a black hole will have its mass increasing 
because of the swallowing of surrounding cosmic matter. 
Our solution represents cases in which there is no accre- 
tion onto the central black hole (as described by = 0) 
and in which the cosmic expansion "wins" over the local 
gravitational attraction. The universe is always expand- 
ing and the density of the cosmic fluid always decreases 
with the cosmic time. The seemingly bizarre behavior of 
the MS mass derives from the fact that it is really a mass 
sum (i.e., the mass of the background fluid is included), 
and it coincides with the Hawking-Hayward quasi-local 
mass [2lj |. In the absence of accretion of cosmic fluid 
onto the black hole, this mass sum decreases in any uni- 
verse in which the cosmological fluid satisfies the weak 
energy condition and its energy density decreases with 
the cosmic expansion. In the presence of phantom matter 
(such as in the late epoch of a quintom-dominated uni- 
verse), the weak energy condition is violated, the phan- 
tom energy becomes more concentrated, instead of being 
diluted, with the cosmic expansion, and the mass sum 
increases. 



Future work will concentrate on studying the general- 
ity of the solutions presented here, and of previous solu- 
tions interpreted as black holes embedded in a cosmolog- 
ical space. To fully support the interpretation of these 
solutions as black holes in a cosmological background, it 
is important that a detailed investigation on the global 
structure of these spacetimes be provided. Until this 
global analysis is available, we cannot yet claim that our 
solutions do represent black holes embedded in expand- 
ing universes, although there is circumstantial evidence. 
This analysis will be carried out in future publications. 
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